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2 ' Abstract. In this paper, we investigate random walks in a family of small- world trees 

^ . having an exponential degree distribution. First, we address a trapping problem, that 

is, a particular case of random walks with an immobile trap located at the initial node. 
We obtain the exact mean trapping time defined as the average of first-passage time 
(FPT) from all nodes to the trap, which scales linearly with the network order N in 
' y ' large networks. Then, we determine analytically the mean sending time, which is the 

mean of the FPTs from the initial node to all other nodes, and show that it grows with 
Q I A^ in the order of A^ In TV. After that, we compute the precise global mean first-passage 

time among all pairs of nodes and find that it also varies in the order of A^ In N in 

the large limit of A^. After obtaining the relevant quantities, we compare them with 

^ I each other and related our results to the eflficiency for information transmission by 

regarding the walker as an information messenger. Finally, we compare our results 
V^ , with those previously reported for other trees with different structural properties (e.g., 

■"sj" ' degree distribution), such as the standard fractal trees and the scale- free small- world 

[~^ ■ trees, and show that the shortest path between a pair of nodes in a tree is responsible 

^D I for the scaling of FPT between the two nodes. 



PACS numbers: 05.40.Fb, 89.75.Hc, 05.60. Cd, 02.10.Ud 
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1. Introduction 

As a powerful mathematical tool that can describe a large number of real natural 
and manmade systems, complex networks have received considerable interest from a 
wide range of scientific communities recently [H Ej |3]. During the last decade, main 
endeavors were devoted to understanding the structural features and dynamics of various 
networks |1]. In particular, treelike networks have attracted renewed attention, because 
the so-called border tree motifs are present in numerous real-life systems and play a 
significant role O [6] . The absence of loops in a treelike network has a drastic influence 
on diverse dynamic processes running on the network, e.g., the voter model [7j and 
naming game [8]. 

Among a plethora of dynamics, random walks in trees have received increasing 
attention in recent years, since the problem is related to a wide range of research fields, 
such as physics [9], biology [10], and cognitive science [H]. One of the most important 
quantities of random walks is the first-passage time (FPT) [12] defined as the expected 
time for a walker starting from a source point to first arrive at a target node [131 fl^ . 
which encodes much information about random-walk dynamics. Thus far, random walks 
in trees with different structures have been intensively studied, including the standard 
fractal trees (e.g., the T-fractal [IS HSl [13 [HI [19] and the Vicsek fractal [201 [21]) and 
scale- free trees [221 ESI El]- These works uncovered how the mean first-passage time 
(MFPT), i.e., the average of FPTs between some given pairs of nodes, scales with the 
network size (number of nodes), and was thus helpful for understanding the impact of 
structural properties on the behavior of MFPT. For example, it was shown that the 
MFPTs between two nodes in different trees behave different. But the main reason for 
the difference remains not well understood. On the other hand, in contrast to the scale- 
free behavior [25], some real networks display an exponential distribution as well [26] . 
Relevant work on random walks on such networks is much less. Particularly, what is the 
main factor affecting the speed of diffusion in general trees is still not well understood. 
A goal of this work is to answer this question, at least partially. 

In this paper, we study a simple random walk [27] on a family of deterministically 
growing small- world trees exhibiting an exponential form of degree distribution [28] . 
We first address a trapping problem, which is a particular random walk with a single 
trap positioned at the initially created node of the networks. We derive analytically 
the mean tapping time (MTT) defined as the average of the FPTs from all nodes to 
the trap, which varies lineally with the network size A^. We then investigate the partial 
mean first-passage time (PMFPT) (i.e., the average of the FPTs from the initial node to 
a randomly selected target node) and the global mean first-passage time (GMFPT) that 
is the average of the FPTs among all pairs of nodes in the networks. Both PMFPT and 
GMFPT are determined through the connection between random walks and electrical 
networks. In contrast to the MTT, both PMFPT and GMFPT are asymptotic to A^ln A^ 
for large networks. We relate our results to the efficiency of information diffusion by 
considering the walker as an information messenger. We also compare our results with 
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Figure 1. (Color online) Construction method of the networks. Each existing node 
(in black) creates a path of m nodes including the mother node itself. 
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Figure 2. (Color online) Illustration for the growing process of a special network 
corresponding to ?7i = 3. 



those found for other trees with different architectures and consequently give possible 
reasons for the behavioral difference of random walks between the considered trees and 
other comparable trees. 



2. Network model and its properties 

We first introduce the network model under consideration, which is built iteratively and 
has a treelike structure [28]. Let Ug [g > 0) be the family of networks after g iterations. 
Initially {g = 0), Uq is a single isolated node without any edge, called the initial node 
below. For g > 1, Ug is obtained from Ug-i by adding a path of m nodes (m is a 
natural number equal to or greater than 2) to each existing node in t/g-i, see Fig. [H By 
construction, it is easy to know that the numbers of nodes and edges in Ug are Ng = m^ 
and Eg = Ng — 1 = m^ — 1, respectively. Figure [2] shows the growing process for a 
particular network for the case of m = 3. Notice that for m = 2, the model reduces to 
the deterministic uniform recursive tree proposed in [29] , which has been extended and 
extensively studied thereafter [301 [SH [321 [33l [M] . 

The networks considered here have some properties observed for many real 
networks. According to the construction process, at any new iteration, the degree 
of every old node increases by 1, independently of the node degrees. Thus, the degree 
distribution of the networks is exponential, instead of a power law [25]. On the other 
hand, the diameter (i.e., the maximum of shortest-path distance between all pair of 
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Figure 3. Another construction of the network. The networks after .g + 1 iterations, 
Ug+i, consist of m rephcas of Ug denoted Ug , t/g , Ug , ■ • •, and Ug , which are 
connected to each other by adding m — 1 edges among the m copies of the initial node. 



nodes) of Ug is {2g — l){m — l) that grows logarithmically with the network size, showing 
that the networks are of small- world [35]. Moreover, some other properties, e.g., the 
adjacency spectrum [28], can also be determined analytically. 

3. Random walks on the networks 

After introducing the family Ug of networks and their properties, we will study the 
discrete random walks [27] performed on Ug. At each time step, the walker (particle) 
jumps uniformly (i.e., with the same probability) from its current position to any of its 
neighboring nodes. One of the most important quantities characterizing such a random 
walk is the FPT [12]. Let Fij{g) denote the FPT for a walker, staring from node i in 
Ug to first arrive at node j. What we are concerned with is how the scalings of FPTs 
behave as the network size increases. 

In the following, we will focus on three cases of random walks. Firstly, we will 
investigate a trapping issue, namely, random walks with a single immobile trap located 
at the initial node, and determine the MFPT to the initial node averaged over all nodes 
in Ug. Then, we will compute the MFPT from the initial node to another node selected 
uniformly from all nodes in Ug. Finally, we will determine the MFPT between all pairs 
of nodes in [/„. 



3.1. MFPT from all other nodes to the initial node 

First, we study a particular trapping problem on Ug, in which the single trap is positioned 
at the initial node. To facilitate computation, we introduce an alterative construction 
method for the networks, which highlights their self-similar architecture, as follows. 
Suppose one has Ug. The next iteration of the network, Ug+i, can be obtained by 
joining m copies of Ug in a way as illustrated by Fig. El 

For convenience of description, we label the initial node in Ug as 1, while the 



duplicates of initial nodes in Uy:i{ri = 2, 3 



{v) 



m. 



are sequentially labeled as 2, 3, 
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m. In addition, let Ti{g) be the FPT, Fi^i{g), also called the trapping time, of node i in 
Ug, which is the expected time for a walker starting from i to first visit the trap node. 
Obviously, for all g > 0, Ti{g) = 0. Let Ttot{g) express the total trapping times for all 
nodes in Ug, i.e., 

Ttot(^) = $^:?^,(^). (1) 

i€Ug 

Then, the MFPT, also called the mean trapping time (MTT) and denoted by {T)g, is 
the average of Ti{g) over all starting nodes distributed uniformly in Ug, given by 

(T), = ^J:U9) = ^. (2) 

Thus, to obtain {T)g, one should first determine Ttot(5')- 

We first compute the FPTs, -Fj+i^j {i = 1, 2, . . . , m — 1) between two arbitrary copies 
of the initial node of Ug^i that are directly connected to each other. Since Ug has a 
treelike structure, according to the result obtained previously in [HI |36l |37] , e.g., Eq. (2) 
in [18], we have 



F,+,,{g) = 2{m - z)m3~' - 1 (3) 

for alH G [1 , rra — 1] , yielding 

Ttot(^) = Ttot(^-i) + $^ Y, Ug) 

m 

= mTUg - 1) + 5^[iV,_iF,,i((7)] . (4) 

4 = 2 

Considering the treelike structure of Ug, we have 

FiM = F,,,_i((7) + F,„i,,„2(^) + ■ ■ ■ + F^,2{g) + F2M . (5) 

Plugging Eqs. ([3]) and ([5]) into Eq. (j4]), we obtain 

Ttotig) = mTtotig - 1) + -(2m2 - 3m + l)m^3-^ - 7^(m - l)m' ■ (6) 

With the initial condition rtot(O) = 0, Eq. (Q is inductively solved, giving 



^tot(5') = m^'^ l:^m - -j -m^ 



^^2\ (g I 

- + -m— — h- 
2 3/ V2 3 



(7) 



Inserting Eq. ([7]) into Eq. ([2]), we obtain a closed- form expression for the MTT on 
network Ug, as 

(r>,^™.(H™-l)-(|.|),„.(|4). (8) 

Recalling Ng = m^, we have g = In Ng/ Inm. Thus, {T)g can be expressed in terms of 
the network size A^^ as 

(r>,^(|-«4)A',-^(m-l)-(|m4). (9) 



Random walks in small-world exponential treelike networks 6 

Therefore, for large networks (i.e., Ng — > oo), 

{T), ~ Ng , (10) 

implying that the MTT (T)^ increases linearly with the network size, independently of 
the degree of the initial node. This linear scaling is in sharp contrast to that of random 
exponential trees [23], in which the MTT depends on the degree of trapping node. 

3.2. MFPT from the initial node to all other nodes 

By definition, Fii{g) denotes the FPT of the walker visiting node i for the first time, 
assuming that the walker started at the initial node in Ug. Let {H)g represent the mean 
value of Fi^i{g) averaged over all target nodes i in network Ug, called the partial mean 
first-passage time (PMFPT). Then, {H)g is given by 

Thus, the problem of finding {H)g is reduced to determining the sum ^^^^ Fi^i{g), 
denoted i^tot(5')- 

Unfortunately, the method for computing Ttot(fi') is not suitable for i^tot(fl')- So, 
we seek for a feasible technique to derive i^tot(fi')- Below, we will apply the hnk 
between effective resistance and the FPTs for random walks [3HIE9] to calculate if tot (5') 
analytically. For this purpose, we replace each edge of Ug by a unit resistor to obtain the 
corresponding resistor networks. To do so, let Ri,j{g) represent the effective resistance 
between two nodes i and j of Ug. Then, we have [381 [39]: Fi,j{.9) + ^jAd) — '^^g^,ji.9)i 
which leads to 

H,ot{9) + Ttot(^) = 2Eg J2 RiAa) ■ (12) 

i&Ug 

Equation (fT2|l shows that if we have the sum Xliec/ -^i,«(fi') °^ ^^^ right-hand side, 
then we can easily obtain Htot{g)- Since, for any tree, the effective resistance Ri,j{g) is 
equal to the geodesic distance dij{g) between i and j, this makes it possible to determine 
the sum ^jg^ RiAq)- To this end, we introduce a new quantity dg, which is the sum 
of shortest distances between the initial node 1 and all other nodes in Ug. By definition, 
we have 

dg = Y.^^M- (13) 

idUg 

Considering the self-similar network structure (see Fig. [3j), we can easily obtain the 
recursion relation 

dg=m dg_, + —-^ '-Ng.l . (U) 

Using di = m{m — l)/2, Eq. (I14p is solved, giving 

dg = ^g{rn - l)m^ . (15) 
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^tot(5') = '^Egdg - Ttotig) 



m 



29 



X 2 1 

g{m-l) -3"^+ 3 



m^ 



m- 1 2 1 

q H — m. 

2^3 3 



and 



{H)g = m^ 



1^ 2 1 

g{m- l)--m + - 



m-1 2 1 

2^3 3 



Equation f ITTj) can be rewritten as a function of the network size Ng, as 

[(m.-^^]-nN. 2 11 r rm, - nin A^. 

{H)9 = N, 



(16) 



(17) 



m-1 lniV„ 2 1 

^ —^ ^ m + - 

mm 3 3 



m-1 lniV„ 2 1 

-^ ^ + -m 

2m?7i 3 3 



Therefore, in the hmit of the large network size Ng, 

{H)g^NglnNg. 



(19) 



3.3. MFPT between all node pairs 

In what follows, we will calculate the MFPT {F)g among all node pairs in Ug, commonly 
called the global mean first-passage time (GMFPT) [ID]. It should be noted that the 
GMFPT has been studied in [31]. Here we will derive an equivalent result using an 
approach different from but relatively easier than that in [34] . 
By definition, {F)g is given by 



F ( ) 1 

\^'9 = AT (AT TT ^ T^T (Aq _ ]_N Z_^ Z_^ ^ijid) : 



Ng{Ng - 1) Ng{Ng - l) 



ij^j i=i 



where the sum 



(20) 



(21) 



denotes the sum of FPTs among all pairs of nodes. Hence, all that is left to find {F)g 
is to determine Ftot(fi')- 

According to the relation between FPTs and the effective resistance, we have 



FU9) = EgJ2J2^^M- 



(22) 



For brevity, we use Dg to denote Xli/j Xljii ^ijid)^ which is the total geodesic distance 
among all pairs of nodes Ug. 

Since Ug can be obtained by the juxtaposition of m copies of Ug-i (i.e., Ug\, Ug\, 
■ ■ ■, and Ug^l) at the edge nodes (replicas of the initial node in f/g-i), Dg can be recast 
as 



Dg = mD„_i + A 



9 1 



(23) 
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where A^ is the sum over all shortest paths whose endpoints are not in the same copy 
off/3-1. 

Denote A"''^ as the sum of all shortest paths with endpoints in f/^_i and f/f_i, 
respectively. According to the value of the distance between two edge nodes in t/°_i 
and tv_i, we can partition the sum of path length A"'^ into m — 1 classes, denoted by 
A"'^(g) with g = l,2, ...,m — 1 being the distance between the two boundary nodes 
in Ug_-^ and Ug_^. That is, A^''^ = X]^ ^'^'^{l)- f^ is easy to see that the number 
of elements in class A"''^(g) is ?7i — g. On the other hand, any two elements belonging 
to A^'^(g) have an identical length of ANg^idg^i + 2q{Ng^iY. Then, the total crossing 
path length A^ can be expressed as 

m—l 

Dg = mDg., + ^(m - q)[ANg_,dg^, + 2g(iV,_i)2] . (24) 

With the initial condition Di = m{m — l){m + l)/3, Eq. ( l24|) is solved to yield 

Dg = -mP [(3^ - 2)m^+^ - (3^ - \)mP + 2m - l] . (25) 

Using the above-obtained results, the expression for (-F)g reads 



^ ^' Ng{Ng-l) N, 



m^ 



X 2 1 

g{m - 1) - -m + - 



3 3 

which can be rewritten in terms of A^^ in the following form: 



g 

2 1 

m - -, (26) 



{F), = iV, 



m-1 lniV„ 2 1 

^ ' ^ m+ 



5" 4- (^^) 



Inm 3 3 

consistent with the result previously obtained in 

Equation f l27|) uncovers the explicit dependence relation of {F)g on the network size 
Ng and the parameter m. In the case of A^^ —t- 00, we have the following expression: 

(F),~A^,lnA^,. (28) 

3.4- Analysis and comparison 

Our results can be related to the efficiency of information transmission. Notice that 
if we consider the walker in the random-walk dynamics as an information messenger, 
then the {T)g measures the efficiency of the initial node (as a receiver) in receiving 
information, while {H)g shows how efficient of the initial node is as a sender to transmit 
information to other nodes, and {F)g is the efficiency of information sending when the 
sender is distributed with equal probability among all nodes. 

The above results, provided in Eqs. (fTOj) and (TT9|) . show evidently that the dominant 
behaviors for (T)g and {H)g are different. The former follows {T)g ~ A^, while the latter 
obeys {H)g ~ A^ InA^, greater than the former. This means that the initial node is 
more efficient in receiving information than sending information. On the other hand. 
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Equation (125]) together with Eqs. f IlOp and (1191) means that although the efficiency of 
the initial node in receiving information is higher than that of the average over other 
nodes, its ability of sending information is similar to the others. Equations ( ITO|) . ( TT9|) 
and (128|) also show that the linear scaling of the efficiency for the initial node measured 
by the MTT is not a representative property of the networks, but the A''^ In A''^ behavior 
for the initial node sending information is so. 

Our obtained results can be compared with those previously reported for other 
treelike networks. Equations ( ITOl) and ( l28l) imply that the position of the trap 
significantly affects the scalings of the MTT for the trapping problem with a single 
trap. This behavior is similar to that in the small- world scale-free trees [2H [29], but is 
in contrast to that of trapping in the T— fractals |T6l [18], [19] and the fractal scale-free 
trees [24], where the MTT does not depend on the trap location. In contrast, as shown 
in Eqs. ( [T9l) and ( 128]) . the initial node has the same dominant scaling of the PMFPT 
as that of the average of PMFPTs over all senders. This equality between PMFPT 
and GMFPT among all node pairs has also been observed for other trees, including the 
T— fractals [TG] [TU [19] and fractal scale- free trees [2lj. However, the scaling of PMFPT 
for different trees may have different behaviors. 

In addition, the leading asymptotic Ng In Ng dependence of GMFPT with the 
network size is also compared to the scalings found from other treelike networks 
with different degree distributions. In the standard fractal trees, such as the T- 
fractals [T6l[T8l[T9] and the Vicsek fractals [21], the GMFPT {F) increases superlinearly 
with the network size iV, which has also been observed from the family of scale-free trees 
with fractality [23]. For star graphs, the GMFPT {F) grows linearly with A^ [21]; while 
for linear chains, {F) scales as a square root of A^ [21]. However, for the class of scale- 
free small- world trees [M] [29], the GMFPT (F) also changes with A^ as (F) ~ A^ln A^, 
which follows the same scaling as that of the exponential trees studied here. 

Finally, combining the present work and the previous studies, it can be seen that 
random walks in trees display rich behaviors in the context of the FPT. At ffist sight, 
degree distribution is perhaps the root responsible for the rich phenomena. However, 
in [23], it was shown that the FPT in fractal and non- fractal scale- free networks 
may exhibit quite disparate scalings, meaning that degree distribution alone cannot 
determine the FPT for random walks on trees. We argue that the FPT on trees is 
determined by the short-path length from the resource node to the target node, while the 
impact of other structural properties is encoded in the short-path length, since the FPT 
Fij from an arbitrary node i to j is actually related to the FPTs of those node pairs for 
two directly connected nodes along the unique short-path direction to the target node, 

i.e., Fij = Fii^+Fi^i^ + Fi^i.^^ \-Fi^i^ + Fi^j provided that i-ii- i2- is ix-iy-j 

is the shortest path from i to j. For example, the FPT on the non-fractal treelike 
scale-free networks [23] and on the exponential networks studied here display similar 
behaviors, since both types of networks are of small world with the average distance 
increasing logarithmically in the network size, in spite of that they have distinct degree 
distributions. As another example, the FPT on some standard fractal trees (e.g. the 
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T- fractals [TSl[Tn] and the Vicsek fractals [21]) displays a superlinear dependence on the 
system size, which is also due to their average distance. 

4. Conclusions 

We have presented a detailed analysis of the simple random walks on a class of 
treelike small-world networks exhibiting an exponential degree distribution. We first 
investigated the trapping problem, focusing on a peculiar case with the trap fixed at 
the initial node, and obtained the exact solution to the MTT, the dominating scaling of 
which varies lineally with the network size A^. We then studied the random walks staring 
from the initial node, and determined analytically the PMFPT from the initial node 
to all other nodes, whose dominant behavior scales with N as A^lnA^. Moreover, we 
determined explicitly the GMFPT among all node pairs and showed that the GMFPT 
also increases with N approximately as N In N . We finally related our results in terms 
of information transmission by regarding the walker as an information messenger, and 
compared them with those previously reported results for other treelike networks with 
disparate topological properties. Our work provides new and useful insight into random- 
walk dynamics running on treelike networks, and could further deepen our understanding 
of random walks on a tree [ID] . 
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